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Are there MACHOs in the Milky Way halo? 
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Abstract. Microlensing searches aim to detect compact halo dark matter via its gravitational lensing effect on stars within the 
Large Magellanic Cloud. These searches have led to the claim that roughly one fifth of the galactic halo dark matter may be in 
the form of compact, solar-mass objects. We analyze this hypothesis by considering the goodness-of-fit of the best-fit halo dark 
matter solutions to the observational data. We show that the distribution of the durations of the observed microlensing events is 
significantly narrower than that expected to result from a standard halo lens population at 90 to 95% confidence, casting doubt 
on the lenses constituting halo dark matter. This conclusion may possibly be avoided if (i) the Milky Way halo is sufficiently 
nonstandard or (ii) a large fraction of the events are due to non-halo populations with event durations coincidentally close to 
those of the putative halo population or (iii) individual event durations have been seriously underestimated due to blending. 
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1. Introduction 

Massive compact halo objects (MACHOs) with mass in the 
range 1O~ 8 M to 10 3 M G can be detected via the temporary 
amplification of background stars which occurs, due to grav- 
itational microlensing, when the MACHO passes close to the 
line of sight to a background star (Paczynski 1986 ). In the early 
1990s several collaborations began monitoring millions of stars 
in the Large and Small Magellanic Clouds (LMC and SMC) Q, 
in order to search for Milky Way (MW) compact halo dark 
matter, and indeed a number of ca ndidate microlensing events 
have been observed (Alcock et al. 2000b; Lasserre et al. 2000; 



The EROS collaboration have five candidate events, result- 
ing from three years of observation of 25 million stars in the 



LMC (Milsztajn & Lasserre 2001 ). Due to the small number of 



events, and poor fit to microlensing of several of these events, 
they use these events to place constraints on the halo fraction as 
a function of lens mass, assuming a standard Maxwellian MW 
halo. In contrast, the MACHO collaboration claim that roughly 
20% of the halo is in compact objects with M ~ O.5M , using 
their 13/17 candidates events which result from the 5.7 years of 
observation of 10.7 million stars (Alcock et al. 2000b I. While 



Milsztajn & Lasserre 2001). The interpretation of these mi- 
crolensing events is a matter of much debate. Though the lenses 
may easily reside within the MW halo and constitute halo dark 
matter, it is also possible that the contribution to the lensing rate 
due to other, so far unknown non-halo po pulatio ns of object s 
has been underestimated (see e.g. Bennett 1998; Zhao 1999). 



not in direct conflict with the EROS exclusion limits these 
MACHO best fit parameters lie just outside the EROS exclu- 
sion limits. These differences could be due to EROS covering a 
larger solid angle than MACHO, such that the lensing rate due 
to lenses in the LMC itself should be smaller, as well as the 
use of less crowded fields by EROS as compared to MACHO, 
simplifying the estimation of the event durations. See Milsztajn 
(2001) for an extended discussion of the current observational 



In particular the lensing rate due to stars, or MACHOs, within 
the LMC itself ( Wu|l994| ; Auborg et al. ^999\ Salati et al ~ 



Evans & Kerins [2000fc Gyuk, Dalai & Grie st [2000|), or a very 
thick Milky Way disk (Gyuk & Gates 2001 ) may be significant. 



1999; 
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situation. 

In this paper we revisit the analysis of the MACHO collab- 
orations candidate microlensing events (Alcock et al. 2000b). 
We contemplate whether the advocated interpretation of the mi- 
crolensing observations in terms of a MACHO halo dark mat- 
ter component with Maxwellian velocity distribution within a 
standard MW halo and with best-fit typical mass 0.6 M e and 
halo dark matter fraction / « 20% actually provides, in abso- 
lute terms, a good fit to the data. 
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2. Candidate events 

The MACHO collaboration apply two sets of selection criteria 
to their data. Criteria A is the most restrictive, only accepting 
events with a single highly significant bump in their lightcurve 
as expected for simple, single lens, events, and results in 13 
events. Criteria B is looser with 17 events passing the selec- 
tion criteria, being designed to also accept low signal to noise 
events and exotic events, such as those where the lens is binary 
and the lightcurve has distinctive caustic features. Marginal 
events which are suspected of being supernovae occurring in 
host galaxies behind the LMC are rejected from criteria A but 
kept in criteria B. Further study of event number 22 []has found 
that the source is extended and has emission lines uncharacter- 
istic of a stellar object and it is hence very unlikely a microlens- 
ing event (Alcock et al. 2001a). This extraordinarily long event 
will henceforth be excluded from the sample, leaving 16 can- 
didate events for criteria B. 

Known stellar populations within the MW disk and bulge, 
as well as within the LMC disk, are expected to account for 
approximately one quarter of the observed events (Alcock et 
al. |2000b| >. Follow-up observations of fields containing each 
of the source stars have been carried out using the Hubble 
Space Telescope Wide Field Planetary Camera 2 (Alcock et 



al. 2001t). In the case of event number 5 (which passes both 



sets of selection criteria) a faint red object has been detected 
close to the source star. A chance super-position is extremely 
unlikely and a combined analysis of the HST and microlensing 
data finds that this object is most likely an M dwarf located in 
the MW disk. 

It is possible to obtain more information about the proper- 
ties of a given lens if either the lens or source is within a binary, 
as the microlensing lightcurve then exhibits additional features. 
The lens responsible for event 9 is a binary, producing distinc- 
tive caustic features in the lightcurve from which the source 
crossing radius and henc e the l ens projected velocit y could be 
estimated (Bennett et al. |l996| ; Alcock et al. |2000a[ ). The low 
projected velocity found suggests that either the lens resides in 
the LMC disk or the source star is also binary. Accurate pho- 
tometry of event number 14 (which passes both sets of selection 
criteria) revealed a periodic modulation in the lightcurve. Fits 
to a rotating binary source have been performed which indi- 
cate that it is most likely that the lens resides in the LMC disc, 
though direct spectral observations of the source are required 
to confirm this (Alcock et al. 2001c). 



While there is clear evidence that event number 22 is not 
due to microlensing and that event 5 is due to a non MW halo 
lens, the locations of the lenses responsible for events 9 and 14 
are not unambiguously known. Furthermore no conclusions on 
the locations of the entire lens population can be drawn from 
the locations of the lenses responsible for events 9 and 14, due 
to the biases involved in observing exotic events (Alcock et 
al. 1200 let Honma |1999|). 



A further complication is the estimation of the event dura- 
tions. In crowded stellar fields it is not possible to resolve in- 
dividual stars and some fraction of the baseline flux may come 



from unlensed stars. This phenomena, known as blending, re- 
sults in an underestimate of the event timescale. The MACHO 
collaboration take this into account by carrying out sophisti- 
cated Monte Carlo simulations to find the average factor by 
which the real timescale is underestimated, and then correcting 



each of the fitted durations by this factor (Alcock et al. 2000b 
Alcock etal. |2001ai. 



3. Likelihood analysis 

The MACHO collaboration defines the likelihood of a given 
model as the product of the Poisson probability of observing 
A^ bs events when expecting N exp events and the probabilities of 
finding the observed durations rj (where j - 1, A^bs) from 
the theoretical duration distribution (Alcock et al. 1996j ; Alcock 



et al. |1997| ): 



£ = exp(-iV exp )n^ / / j , 

where the expected number of events is given by 



T 

Jo 



A'oxp = E | ^e(f)df, 



and fij by 



Mj = E e(fj) 



dF(fj) 
df 



(1) 



(2) 



(3) 



2 We use the M ACHO collaboration's event numbering nomencla- 
ture (Alcock et al. pOOObl) throughout. 



Here E = 6.12 x 10 7 star-years is the exposure, e(f) is the detec- 
tion efficiency, t is the estimated event duration statistically cor- 
rected for blending and dT/df is the differential event rate. The 
likelihood approach can be used to compare how well different 
sets of parameter values within a given theoretical model fit the 
data. In the present context these parameters are the MACHO 
mass and halo fraction and the model is a standard halo pop- 
ulation of compact dark matter. However, the likelihood ap- 
proach does not tell one, in absolute terms, how good a fit re- 
ally is. In particular, it is possible that though "best-fit" model 
parameters have been found, the model provides a generally 
bad description of the observational data. In order to illustrate 
this we compare dF/df of the best fit standard halo model with 
delta MACHO mass function (S-DMF) with that of a best fit 
(i.e. maximum likelihood) Gaussian differential event rate with 
variable amplitude, center and width. We note that the Gaussian 
has no physical motivation and is intended to simply provide 
some arbitrary reference to the S-DMF model. 

The best fit Gaussian and S-DMF differential event rates 
which arise in a likelihood analysis assuming that all 13 events 
found using criteria A are due to MW halo lenses (i.e. neglect- 
ing the event rate due to background events) are plotted in Fig. 
1 multiplied by the detection efficiency. The Gaussian clearly 
gives a better fit to the observed event distribution than the S- 
DMF model, and its likelihood is greater by a factor of 32. 
We have also carried out this comparison for all possible 10 
event sub-sets (i.e. subtracted every possible combination of 3 
background events) and in every case the Gaussian has a larger 
maximum likelihood (by a factor ranging from 5 to 4500). 

3 Note that the Gaussian model has one more fitting parameter than 
the S-DMF model. 
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Fig. 1. The best fit delta function plus standard halo model 
(solid line) and Gaussian (dotted) differential event rates for 
the 13 events resulting from criteria A. Crosses indicate the 
observed event durations. An analogous plot for the 16 events 
resulting from criteria B appears very similar. 

4. Event duration distribution 

Fig. 1 illustrates the reason why the artificial Gaussian event 
duration distribution provides a better fit to the data than the 
event duration distribution produced by a MACHO population 
within a standard halo model; namely, whereas the S-DMF pre- 
dicts the occurrence of long duration events none are observed. 
To quantify this discrepancy we have performed a Monte Carlo 
simulation which compares the observed width of the distribu- 
tion 



Width 



Max ( j=hN;j(tj) - Min u= i^ fa} (^) 



Avg,, = 



(tj) 



(4) 



to the Width obtained by statistically drawing a number N 1 ™ 1 " 
of events from a standard halo event duration distribution (mul- 
tiplied by the relevant detection efficiency and taking into ac- 
count the transverse velocity of the line-of-sight). We note that 
the above statistic is designed to be independent of the best-fit 
MACHO halo fraction and mass (] and is thus a parameteriza- 
tion of the observations which is essentially orthogonal to the 
best-fit parameters determined by the MACHO collaboration^]. 

In our Monte Carlo analysis we find the fraction of gener- 
ated event samples that yield a Width which is smaller than that 
observed. To account for background (i.e. non-halo) lenses we 
employ the estimates for the expected nu mber of background 
events < N B c > found in Alcock et al. (|2000b[). We do not 



distinguish between the different background populations, but 
simply subtract a number N* BC from N* b to obtain A 7 *?* , the 
number of events which are expected to be due to MACHOs 
and are to be generated from the theoretical event duration dis- 
tribution. We do this using two different methods: subtracting 

4 There remains a residual dependence of the Width on the typi- 
cal MACHO mass, however, due to the duration dependence of the 
detection efficiencies. 

5 This approach is similar to that of de Rujula et al. (1995) who 



Table 1. Probability of various subsamples of the observed can- 
didate microlensing events from the 5.7 year candidate events 
being drawn from an event duration distribution function (i.e. 
dT/df) as predicted by the standard MW halo model. Here the 
columns indicate the number of candidate events (N„i, s = 13 
for efficiency A, and N bs = 17 for efficiency B), whether a 
delta- (DMF) or Gaussian- (GMF) MACHO mass function has 
been assumed, the numbers of the events which have been ex- 
cluded from the sample by hand (correspondingly reducing the 
number of observed events from N bs to N* b ), the assumed av- 
erage number < Nbc >* of events due to non-halo background 
(BG) lens populations, as well as the aforementioned probabil- 
ities, computed either by subtracting a fixed number < Nbc >* 
from N*. to obtain the number of halo events N llc ! l " to be drawn 

obs obs 

from the theoretical event duration distribution, or by finding 
the number iV* of events to be subtracted from N*. , from a 

BO obs 7 

Poisson distribution with average < Nbc >*■ In the case of non- 
integer < Nbc >*> this number is rounded up to the next integer, 
when fixed number background subtraction is performed. Note 
that the exclusion of event number 22 from the sample by hand 
is not assumed to reduce < Nbc > from 3.9 to 2.9, as it is a su- 
pernova rather than microlensing due to a background lensing 
population. 

Table 1 



Model/ 


mass 


events 


# BG 


Prob. 


Prob. 


Sample 


funct. 


excluded 


events 


(fixed) 


(Poisson) 


S/13 


DMF 


none 


3.0 


6.9% 


8.9% 


S/13 


GMF 


none 


3.0 


3.5% 


5.2% 


S/13 


DMF 


5 


2.0 


6.9% 


9.1% 


S/13 


GMF 


5 


2.0 


3.5% 


5.2% 


S/13 


DMF 


5, 14 


1.0 


8.1% 


8.7% 


S/13 


GMF 


5, 14 


1.0 


4.2% 


4.7% 


S/17 


DMF 


22 


3.9 


6.4% 


7.7% 


S/17 


GMF 


22 


3.9 


3.1% 


4.2% 


S/17 


DMF 


5, 9, 14, 22 


0.9 


4.7% 


4.8% 


S/17 


GMF 


5, 9, 14, 22 


0.9 


2.2% 


2.3% 



a fixed number A^ G =< Nbc >*> an d subtracting N* BC1 as drawn 
from a Poisson distribution with average < Nbc >*■ [J We note 
here that while the background events could potentially account 
for the longest or shortest events in the sample, we nevertheless 
do not change the observed Width on subtracting background 
events, such that our estimated probabilities, in a conservative 
spirit, should generally reflect an overestimate. 

The results of this procedure, in particular the fraction 
of simulations which have smaller Width than observed, are 
shown in Table 1 for a variety of different subsamples. This 
fraction is generally between 5% and 10%, such that the under- 
lying S-DMF model serving as an explanation for the observa- 
tions is ruled out at the 90 to 95% confidence level. Note that 
this conclusion does not vary strongly if the events which are 
suspected of being due to background populations are excluded 
by hand. This would not be the case if event number 22, which 
has a duration of 297.8 days, was a genuine microlensing event. 



showed that the duration dispersion could in principle be used to sep- 
arate background events. 



6 In some of our samples < N BG > as given in Alcock et al. ( 2000b| ) 
is adjusted to < N B o >* to account for events which have been ex- 
cluded from the sample by hand, cf . Table 1 . 
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The discrepancy between the observations and the theoret- 
ical expectations even worsens when the assumption of a delta 
function MACHO mass function (DMF) is dropped. To illus- 
trate the effects of a MACHO mass function with finite width 
we take the number of MACHOs with mass within a logarith- 
mic interval to be 



dP 



1 



dln(M/M*) 
with 



V2^ 



i / uHw) 



(T 



cr 



cr= Vln2/2, 



(5) 



(6) 



such that MACHOs with either twice or one-half the average 
mass < M >- exp(<x 2 /2)M* (which is assumed to take the 
best-fit value of the DMF) are a factor of two less common. 
These models are indicated in the table by GMF. 

It may be noted that the probabilities are generally larger 
when background event subtraction is performed by employing 
Poisson statistics, rather than by subtracting a fixed number of 
events. This is due to the probabilities in the Poisson case being 
dominated by samples which have N* BC larger than < Nbg >*> 
and have correspondingly smaller N^'°, making a comparison 
of theoretical expectation and observed distribution more fa- 
vorable. Nevertheless, such reasoning makes the implicit and 
nontrivial assumption that all background events have event du- 
rations falling within the range of those of a putative MACHO 
population. 

One may ask, making this assumption about coincidence 
in event durations between different populations, by how much 
one has to increase the number of background events, in order 
to achieve a "decent" probability (say 20%) of the observed dis- 
tribution indeed being drawn from the distribution expected for 
the standard halo model. Fig. 2 shows the probability for differ- 
ent samples, for the delta mass function (DMF) and that given 
by Eq. 5 (GMF), indicating that an increase of N* BG by about 
a factor of two would be necessary in order for this to be the 
case. This suggests, on statistical grounds, that the contribution 
of the background lens populations to the observed events has 
been underestimated, and consequently the MACHO halo mass 
fraction overestimated. Nevertheless, due to the small number 
of observed events, a significant contribution to the observed 
events by a population of MACHOs within a standard MW halo 
may not be ruled out. This is interesting in light of the recent 
detection of a population of old very cool white dwa rfs wh ich 
could comprise a few per-cent (Opp enheimer et al. 2001 ) of 
the MW dark halo, see Richer ( 2001 ) for a critical discussion 
of these observations. 



5. Conclusions 

In this paper, we have shown that the most straightforward in- 
terpretation of the existing microlensing observations of stars 
within the LMC, assuming the majority of candidate events to 
be due to compact lenses within the MW halo, may actually 
be in conflict with the observational data itself. The MACHO 
collaboration uses a likelihood analysis to derive two "best- 
fit" parameters from their microlensing observations: the typ- 
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Number of background events 
Fig. 2. Probability of obtaining an event duration distribution 
for a standard halo model as narrow as that observed in the 
5.7 year sample as a function of the number of events as- 
sumed to be due to background population(s). Here the back- 
ground events are implicitly assumed to neither contribute the 
longest nor shortest event, but to fall, by pure coincidence, 
within the range of the event durations of the putative MACHO 
population. Background subtraction is performed by subtract- 
ing a fixed (non-Poissonian) number of events from the num- 
ber of observed candidate events. The lines show results for 
model/sample S/13 - DMF (solid), S/13 - GMF (dotted), S/17 - 
DMF (dashed), and S/17 - GMF (dashed-dotted), respectively, 
where in the latter two samples event number 22 has been ex- 
cluded. 



ical MACHO mass and MACHO halo fraction. A naive com- 
parison of the theoretical event rate distribution produced by 
these "best-fits" assuming a standard halo model, with an ad- 
hoc Gaussian event rate distribution indicates that the spread in 
the durations of the observed events seems smaller than would 
be expected for a standard halo of compact dark matter. In or- 
der to quantify this suspicion we have formulated a statistic, 
the Width (essentially the spread in the observed event dura- 
tions relative to the average event duration), that is virtually 
independent of the MACHO best-fit parameter values. 

We have carried out Monte Carlo simulations to calcu- 
late the expected distribution of the Width, under the assump- 
tion of a standard halo MACHO population, utilizing two 
different background subtraction methods (fixed number and 
Poissonian), for both detection criteria, for delta function and 
finite width MACHO mass functions and excluded various sets 
of events suspected of being due to background populations. 
The probability of observing a value of the Width as small as 
observed falls, in all cases, below 10%, with values often be- 
low 5%. In other words the underlying model, in terms of the 
expected background event rate and a standard halo MACHO 
population, is excluded at between the 90% and 95% confi- 
dence level. 

The discrepancy between the observational data and its halo 
dark matter interpretation may be alleviated if the number of 
events which are due to non-halo populations is significantly 
larger than expected, at the expense of correspondingly reduc- 
ing the MACHO halo fraction. Such a reasoning also implic- 
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itly assumes the non-halo population to produce microlensing 
events with event durations being by chance very similar to 
those of the halo population. Possibly more natural would be 
an explanation in terms of "all background", though the na- 
ture of this background, and the reason for it producing such 
a narrow distribution in durations, remains to be established. 
Alternatively, the discrepancy may indicate that the standard 
halo model is a poor approximation to the actual MACHO dis- 
tribution (Widrow & Dubinski 1998 ). For instance, there could 



be a clump with small intrinsic velocity dispersion crossing 
the line of sight to the LMC, which, nevertheless, would be 
required to be surprisingly massive. Another possible resolu- 
tion to the problem could be that the statistical correction to 
the event durations, due to blending, leads to a reduction in 
the spread of the durations i.e. that the underlying durations of 
some observed events are significantly longer than their esti- 
mated values. If this were the case, however, then conclusions 
based on these estimated timescales, not only the analysis in 
this paper, but also the derivation of exclusion li mits on the 
halo fraction in massive black holes (Alcock et al. 2001a ) due 
to the purported absence of long-duration events, would seem 
unreliable. 

Further clarification of the nature of the lenses may hope- 
fully result from an enlarged sample of microlensing candi- 
dates, as expected from the upcoming analysis of the MACHO 
8 year data and EROS 5 year. Evidently a longer survey du- 
ration will increase the number of candidate events, reducing 
the ambiguity associated with Poisson statistics, as well as in- 
creasing the efficiency for detecting longer duration events, 
thus making the width of the distribution a significant discrim- 
inator between models. Nevertheless, a signifcant reduction of 
Poisson ambiguity may only result by an increase of the num- 
ber of events by an order of magnitude, such as anticipated in 
a next generation SUPERMACHO survey (Stubbs |1998[ ). This 
later survey is expected to also reduce uncertainties in event 
duration corrections due to blending as it is planned to employ 
the technique of difference image analysis. 
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